library(tidyverse)

mpileup <- read_delim('Archive/covid19/data/vindi_ldblock.mpileup.txt.gz', col_names = c('chr','position','ref','depth','reads','quality'))

mpileup |>
  filter(position == g2r_position) |>
  mutate(ref_reads = str_count(reads, '\\.|,'),
         alt_reads = depth - ref_reads)
  
